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To understand the roles of infall and protostellar evolution 
on the envelopes of massive young stellar objects (YSOs). The 
chemical evolution of gas and dust is traced, including infall 
and realistic source evolution. The temperatures are determined 
self-consistently. Both ad/desorption of ices using recent labo- 
ratory temperature-programmed-desorption measurements are 
included. The observed water abundance jump near 100 K 
is reproduced by an evaporation front which moves outward 
as the luminosity increases. Ion-molecule reactions produce 
water below 100 K. The age of the source is constrained to 
t ~ 8 + 4x1 4 yrs since YSO formation. It is shown that the 
chemical age-dating of hot cores at ~ few x 10 3 - 10 4 yr and 
the disappearance of hot cores on a timescale of ~ 10 5 yr is a 
natural consequence of infall in a dynamic envelope and pro- 
[ tostellar evolution. Dynamical structures of ~ 350AU such as 
£N| disks should contain most of the complex second generation 
species. The assumed order of desorption kinetics does not af- 
C3 • feet these results. 

Key words, stars: formation - stars: individual: AFGL 2591 - 
| ISM: molecules - molecular processes 



(N 
> 

m 
m 
in 
o 

O 

6 

in 
c3 



Astronomy & Astrophysics manuscript no. dotyetal 


February 4, 2008 


(DOI: will be inserted by hand later) 





Astrochemical confirmation of the rapid evolution of massive 
YSOs and explanation for the inferred ages of hot cores 

S. D. Doty, 1 E. F. van Dishoeck, 2 and J. C. Tan 3 

1 Department of Physics and Astronomy, Denison University, Granville, OH, 43023, USA 

2 Sterrewacht Leiden, PO Box 9513, 2300 RA Leiden, The Netherlands 

3 Department of Astronomy, University of Florida, P. O. Box 112055, Gainesville, FL, 32611, USA 



Received ; accepted 
Abstract. 

1. Introduction 

The physical and chemical structure of massive young stel- 
lar objects (YSOs) has been an area of increasing study (e.g. 
Beuther et al. 2006 1 While a wealth of new data and models 
have become available, many basic open questions remain. For 
example, the processes and timescales involved in the forma- 
tion and evolution of the high-mass protostellar objects and en- 
suing "hot core" phase are not on firm theoretical ground. 

Recently, McKee & Tan ( 2003 1 introduced the first realistic 
model for the evolution of massive YSOs, including a time- 
dependent accretion rate. By combining this with a model for 
the evolution of the protostar onto the ZAMS, they predicted 
the time-dependent total (source + accretion) luminosity for 
the lifetime of the massive YSO. 

A test of the evolution is expected to be desorption of ices 
from grain mantles, as the evolving source heats the surround- 
ing envelope. Laboratory data using temperature programmed 
desorption (TPD) have been recently reported (see e.g., Fraser 
et al. l200H . These data allow for a detailed analysis of the des- 
orption of realistic ice mantles. 

Models have included some, but not all, of these effects. 
While TPD data were included by Viti et al. J2004t and Lintott 
et al. (2005), these were single-point models and/or used ap- 
proximate source evolution. Source evolution was considered 
by Rodgers & Charnley d2003l and Lee et al. (2004i, but us- 
ing low-mass YSO models. It is common to assume a static 
envelope, constant luminosity, and standard first-order desorp- 
tion (e.g., Casel li et al. fl993l Nomura & Millar I20041 Viti & 
Williams ( 1999 1 and others have used multi-point models. 

In this letter we describe the first model for the chemical 
evolution in a collapsing massive YSO, including the evolution 
of the central source, infall, and ad/desorption of ices from the 
grains. The models are used as a clock for studying the evo- 
lution of the central source. Here we concentrate on the distri- 



Send offprint requests to: S. D. Doty, e-mail: doty@denison.edu 



bution of H2O, since it is not only a dominant ice component 
and important thermal regulator, but also a target of the upcom- 
ing Herschel Space Observatory. We apply the model to the 
AFGL 2591 envelope as it is one of the best observed sources 
for chemistry in general, and water in particular. 

2. Source and Model 

AFGL 2591 has been observed in water lines over a range 
of excitation conditions with both S WAS (Snell et al. |2000 1 
and ISO (Boonman & van Dishoeck 2003 1. Boonman et al. 
(2003) used radiative transfer models to constrain the wa- 
ter distribution, with three conclusions: (1) water ice evapo- 
rates at r ~ 10 I6 ~ 165 cm [T ~ 100 K in the static model]; 
(2) x(H 2 0)r>ioo = n(H 2 0)/n(H 2 )r>ioo = 2 x 10" 4 ; and (3) 
jc(H 2 0):t<ioo < 10~ 8 . This distribution has been confirmed by 
recent H^O observations by van der Tak et al. ( 2006 1. 

Van der Tak et al. (2000 1 constructed spherically- 
symmetric models of AFGL 2591 and found a best fit density 
power law n(r) oc r~ m where m ~ 1 . Since the source evolu- 
tion models of McKee & Tan ( 2003 ) break down for m — 1 
exactly as the pressure gradient vanishes in this limit, we adopt 
m = 1.1. We take a final stellar mass of m t j = 20 M . All 
other inputs are taken from McKee & Tan (2003) for cores. 
Consistent with m ~ 1, we adopt a logotropic collapse velocity 
field (McLaughlin & Pudritz 11997b . While chemical and ra- 
diative transfer modeling suggest m = 1 , we also consider an 
m =1.5 polytrope. 

Figure^shows the luminosity evolution, L(t). The results 
are almost independant of m t j. Given n(r, t) and L(t) we solve 
for T(r, t). Figure|5]shows T^ ust (r, f). 

The chemical model is based upon the best-fit model of 
Doty et al. (2002) for AFGL 2591, using the same gas-phase 
chemical network, cosmic -ray ionization rate, and initial abun- 
dances in the cold cloud exterior. In particular, since L ~ at 
t — 0, the input mantle composition places all H2O and CO 
on the grain surfaces intially. The CO is quickly liberated - 
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Fig. 1. Evolution of L(t). The curves within each line type cor- 
respond to final stellar masses of 40,20, 12, and 8 M (top to 
bottom). 

by 3 x 10 3 yr in the H2O evaporation zone, and by 3 x 10 4 yr 
throughout the envelope. 

We have made two important modifications. First, we al- 
low the material to infall, solving the chemical evolution in 
a Lagrangian frame, while simultaneously interpolating on 
T(r, f). Second we include ad/desorption. For water, we incor- 
porate both the zeroth-order kinetics of Fraser et al. ( 2001 1, and 
the first-order kinetics of Sandford & Allamandola (19883, m 
turn. For other species, we follow the usual first-order kinetics. 

3. Results 

The water abundance after including infall and source evolu- 
tion for m = 1.1 is shown in Fig. [3] The abundance constraints 
are noted by the shaded regions, while the region of the ob- 
served evaporation step is noted by the double arrows. 

These results can be understood given the fact that for 
Tdust > 100 K, water ice evaporates quickly. As a result, the 
Tdust > 100 K region lies behind an outward propagating evap- 
oration front at which the water ice evaporates nearly instanta- 
neously from the grain surface. The evaporation front stops at 
r ~ 3 x 10 16 cm and t ~ 2 x 10 5 yr, as m„ — > m m j. For the 
m — 1 .5 poly trope, this occurs at t ~ 9 x 10 4 yr. 

The water abundance for 7d us t < 100 K can be understood 
in terms of ion-molecule chemistry. Here CO is dissociated into 
O. Proton-transfer processes this into H30 + , which then dis- 
sociatively recombines into H2O. The net result is that CO is 
gradually transformed into water over time in the cool exterior. 

This combination of effects brackets the observational con- 
straints shown in Fig. [5] For t < 3 x 10 4 yr, the luminosity 
has not increased to where the evaporation front can match 
the observed size of the region of enhanced gas-phase water 
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Fig. 2. Temperature distribution as a function of position and 
time for m = 1.1, The curves are labeled by the time in yr. 

of r ~ 10 16 ~ 16 5 cm. For t > 10 5 yr, the abundance in the 
outer envelope is too high primarily due to ion-molecule pro- 
duction of water for T < 100 K. Together this means that a 
time t — 3 - 10 x 10 4 yr since YSO formation is sufficient to 
liberate water, without warming the entire envelope or produc- 
ing significant ion-molecule chemistry. Results for an m — 1 .5 
polytrope are nearly identical. 

Lower stellar masses yield commensurately longer times 
to reach a given desorption radius, due to their longer Kelvin- 
Helmholz timescales. Fig. [2 shows that the observational dif- 
ferences are small: L(t) varies by < 0.5 dex for 8 < m % AM Q ) < 
40. Since desorption occurs around hook ~ 4 x 10 14 cm 
(L/3L )^ 2 , a stellar mass 5 times smaller leads to a desorption 
radius at any time that is only marginally (< 1.8x) smaller. 

To understand other implications, in Fig.0]we plot the infall 
trajectories of parcels. Grain mantles quickly evaporate (T > 
Tdes) in the shaded region The lines crossing the trajectories 
signify where dynamics becomes important (Td yn ~ ?)■ 

We first note that for parcels originating at r < 10 16 cm, the 
effects of source evolution dominate collapse, since T reaches 
T^s before dynamics become important. On the other hand, 
parcels originating atr>3-6xl0 16 cm are dynamic in the 
mantle evaporation region, due to the fast evolution of the cen- 
tral source. In all cases the gas spends little time above ra es . 

To see this more clearly, in Fig.|5]we plot the time spent 
by parcels above r<j es as a function of their initial position. 
Qualitatively, the time spent in the warm gas can be understood 
in three regimes. (I) At early times and small distances, the en- 
velope is ~ static and evolution is dominated by the source. 
Parcels that begin further out spend more time at high temper- 
atures due to their longer infall distance. (II) An intermediate 
regime exists where both dynamics and source evolution oc- 
cur. The shorter times at high temperatures due to the faster 
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Fig. 3. Water abundance as a function of position in AFGL 
2591. The different line types correspond to different times 
since the formation of the YSO. The shaded areas are the ob- 
servational constraints. 



Fig. 4. Infall trajectories of parcels (bold), crossed by lines not- 
ing where dynamics become important (thin). The shaded re- 
gions signify desorption (dark = logotrope, light = polytrope). 



infall speeds are somewhat balanced by the outward progres- 
sion of the evaporation front. Finally, (III) cold parcels from 
far out warm as they fall into the completely evolved inner en- 
velope. In this case, parcels that start progressively further out 
hit the now stationary evaporation front at progressively higher 
speeds, leading to less time at high temperatures. 

Quantitatively, Fig. [5] shows that parcels spend ~ few x 
10 3 — 10 4 yr as "warm" gas, independent of their starting point. 
Parcels then enter a region of < 200 AU (< 0.2" at 1 kpc). 
Such regions are severely beam-diluted and thus nearly invis- 
ible. As a result, chemical clocks based upon warm gas-phase 
("hot-core") species - especially observed with single-dish in- 
struments - should show ages of ~ few x 10 3 - 10 4 yr. These 
ages do not sample the actual age of the source, but instead the 
residence time in the extended warm gas. The combination of 
infall, evaporation, and source evolution naturally explains the 
observed consistency of chemically inferred ages of hot cores 
between sources. 

Since the outer envelope remains "cold" for ~ 10 5 yr, we 
predict a discrepancy between the ages inferred from hot core 
molecules and from cold gas-phase species. This is consistent 
with the chemically-inferred ages of ~ 4 - 40 x 10 3 yr for hot 
cores, and > 10 5 yr for their cold surrounding envelopes (e.g. 
Hatchell et al. 119951 Millar et al. lT997l . 

Taken together, we propose the following picture of high- 
mass envelope evolution: for t < 10 4 yr, the majority of the 
observable envelope is essentially static, and a cold gas-phase 
chemistry is seen. As the YSO evolves and warms the enve- 
lope, a warm gas-phase "hot-core" chemistry is produced. For 
t < 1-8x1 4 yr (depending upon poly tropic index), a static hot 
core of size r ~ 10 16 cm is produced. After this time, the "hot" 



material is dynamic, while much of the cold material remains ~ 
static. The ice-evaporated hot core chemistry has a characteris- 
tic inferred age of ~ few x 10 3 - 10 4 yr, independent of the ac- 
tual age of the YSO. Finally, once the source reaches the main- 
sequence, it is hot and luminous enough to ionize the envelope, 
leading to the disappearance of the hot core after t ~ 10 5 yr. 

In the scenario above, the infalling material accretes onto 
a central structure smaller than a few hundred AU. This 
may be an accretion disk that is diluted in large single- 
dish observations. Indeed, the disk scale is given by = 
350 AU O6/O.OO7)(m H ,/2OM ) 1/2 , where f3 is the ratio of ro- 
tational to gravitational energy. For high-mass star-formation, 
/3 ~ 0.007 (Pirogov et al. 2003 1. This size is consistent with a 
~ few hundred AU disk inferred by van der Tak et al. d2006> 
toward AFGL 2591. The column density of this region should 
be > 10 22 cirT 2 , suggesting densities > 10 7 cm -3 and tempera- 
tures above 1 00 K. With lifetimes of 1 5 - 1 6 yr (Cesaroni et al. 
2006 1, this is the regime where warm, complex, second gener- 
ation chemistry can take place. As a result, future observatories 
such as eSMA, ALMA, and the EVLA, which can probe scales 
of < 0.2" will be exquisite tools for uncovering and age-dating 
the complex circumstellar chemistry. 

Lastly, the first-order desorption kinetics of Sandford & 
Allamandola ( 1988 1 yield similar results to the zeroth order ki- 
netics from TPD data. The abundance distributions have the 
same shape, but shift inward by ~ 0.1 dex in position due to 
their lower desorption energies. While the order of desorption 
kinetics is not significant here, it may be important elsewhere. 
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Fig. 5. Time spent above 7d es in the extended envelope. 



4. Conclusions 

We have constructed the first model for the chemical evolution 
in a collapsing massive YSO, including the realistic evolution 
of the central source as well as the ad/desorption of ices from 
grain mantles as the grains infall from the cool exterior into 
the warming interior. This approach allows for a more real- 
istic evaporation of ices from the mantles, as opposed to the 
more parametric assumptions of previous models. We find that 
this approach naturally reproduces and explains the parametric 
step-function water distribution inferred observationally. 

The water abundance jump is naturally explained as the 
propagation of a thermal desorption front outward from the 
YSO as it evolves. The abundance in the cool outer envelope 
is explained by ion-molecule reactions which convert CO into 
H2O. The time since protostar formation is bracketed by these 
two processes, yielding an age constraint off ~ 8+4x1 4 
yrs since formation of the YSO. Due to the weak dependence 
of T(r) on L(f), this timescale should be similar (within 0.5 - 1 
dex) for a range of stellar masses. 

The combination of the warming envelope and material in- 
fall leads to small residence times of material in the warm gas, 
~ few x 10 3 - 10 4 yr, after which the material is beam-diluted 
in a small structure such as a disk. This naturally explains both 
the small variance in the ages of hot cores inferred from chem- 
ical clocks, as well as the discrepant inferred ages between hot 
cores and their surrounding halos. Hot cores should disappear 
in ~ 10 5 yr, once the YSO reaches the main sequence and ion- 
izes the envelope. These results suggest that massive star disks 
should contain most of the complex species on scales of < few 
hundred AU, which will be observable by the next generation 
of submm interferometers. 



Finally, the order of desorption kinetics assumed (zeroth 
vs. first-order) are not significant in comparison to source evo- 
lution or observational uncertainties. 
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